Multiscaling in the YX model of networks 
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We investigate a Hamiltonian model of networks. The model is a mirror formulation of the XY model (hence the 
name) — instead letting the XY spins vary, keeping the coupling topology static, we keep the spins conserved 
and sample different underlying networks. Our numerical simulations show complex scaling behaviors, but no 
finite-temperature critical behavior. The ground state and low-order excitations for sparse, finite graphs is a 
fragmented set of isolated network clusters. Configurations of higher energy are typically more connected. The 
connected networks of lowest energy are stretched out giving the network large average distances. For the finite 
sizes we investigate there are three regions — a low-energy regime of fragmented networks, and intermediate 
regime of stretched-out networks, and a high-energy regime of compact, disordered topologies. Scaling up the 
system size, the borders between these regimes approach zero temperature algebraically, but diff'erent network 
structural quantities approach their T = 0- values with diff'erent exponents. 



I. INTRODUCTION 

The XY model is one of the classical and most versatile 
spin models of statistical mechanics. It has been used in many 
physical contexts, from superconductors ([25]) to lattice gauge 
theory of quantum chromodynamics ( Il8| ), and further on ap- 
plied to emergent phenomena in systems like bird flocks (l26l) 
and parallel discrete-time simulations ( 15 ). In general, the XY 
model describes a system of N pairwise interacting units. Let 
the interaction be represented by a graph G = (V,E) where 
V is the set of units, or vertices, and E is the set of edges 
(pairs of coupled units) Each unit i e V is characterized by a 
real number 6i called the spin of /. The probability of a cer- 
tain combination G and {^/j/ey is proportional to exp (-H/T) 
where the parameter T is called temperature (parameterizing 
the disorder of the system), and H is the Hamiltonian 

H= - Yj ^os(6i-6j) (1) 

(i,j)eE 

— a function that in physical systems represents the energy of 
a configuration. (Note that, in our formalism, temperature and 
energy are dimensionless.) From the symmetry of the cosine 
function in Eq. ([T]), we see that the values of Oi only matters 
modulo 27C. For this reason is Oi commonly restricted to the 
interval (0, In]. 

Traditionally, studies of the XY model take G as a lattice 
graph — a graph that can be drawn as a lattice (a discrete sub- 
group spanning the vector space R^, where d is the dimension, 
a natural number) — and considered fixed, while {6i}iEV is the 
object of study. (The most important result is perhaps that 
for a two-dimensional lattice structure the XY model under- 
goes a peculiar phase transition, the Kosterlitz-Thouless tran- 
sition ([16 ) between a disordered phase and a phase with alge- 
braic spin correlations.) Sometimes random graph ensembles 
are used ([T3l [TTl [27]) . sometimes regular topologies that are 
not mathematical lattices ( 3 ), but usually G is static (the only 
exception we are aware of is Ref. (J ), where the XY model is 
simulated on a graph being rewired without any bias). In this 
paper we turn the situation around — we fix {^/j/ev and let G 
be only restricted by N, the number of edges M and that the 
graph is simple (i.e., that there are no multiple edges or self- 
edges). The spins are, like common practice when initializing 



the XY model for Monte Carlo simulations, drawn with uni- 
form probability. Our model is a conceptual mirror image of 
the XY model — hence we call it the YX model. 

The YX model defines an ensemble of graphs and is thus 
more related to the recent works on models of complex net- 
works (f^ ^). Authors have studied general classes 
of graphs sampled with a Boltzmann-like probability — so 
called exponential random graphs (24), Markov graphs (lZj), p- 
star models ( 2 ) or statistical-mechanics models ( [22]) depend- 
ing on the background of the author. These models are used 
for sampling networks with some prescribed structures; they 
are thus not microscopic, or mechanistic, models and usually 
not analyzed in terms of emergent properties in the N ^ oo 
("thermodynamic") limit. The YX model is also related to hid- 
den variable models ([H [T4[ [19]) where some variable (in this 
case 6i) assigned to the vertices is aff'ecting their position in 
the network. 

The main reason of investigating the YX model is somewhat 
academic — what happens when we update the connections, 
rather than the spins, in a simple statistical-physics model? A 
token that the YX model has a rich and interesting behavior is 
that the ground state is frustrated (i.e., the terms, - cos (Oi-Oj) 
in the sum of Eq. ([T]) is not minimal (-1) for all pairs (/, j) e E. 
If we make a similar modified Ising, model or ^-state Potts 
model, if the number of edges is small enough (smaller than 
the number of edges in a graph where all nodes of the same 
spin have an edge between them), the system is not frustrated. 
Such models will most likely have a second-order phase tran- 
sition between segregated and mixed states and might be in- 
teresting in their own rights. In the YX model it is hard to 
guess the outcome, even the ground state, a priori. 

In the rest of the paper we will present the simulation 
scheme in more detail and analyze the size scaling of network 
structural quantities. 



II. SIMULATIONS 

We simulate the YX model using Metropolis Monte Carlo 
sampling (21 ). For one update step, we choose three distinct 
vertices /, j and k such that (/, j) e E but (/, k) ^ E at random. 
Let E' be E with (/, 7) replaced by (/, k), then we accept the 
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(a) -199.52 (b) H^-195M (c) -192.05 




FIG. 1 Three low-energy configurations of the YX model with 100 
vertices. 



change (replacing E by E') if H(E) > H{E')\ or, if H{E) < 
H(E'), with a probability proportional to 

expl 1. (2) 

Let, furthermore, N update steps comprise one sweep. 

The configuration space of the YX model is, as will be dis- 
cussed later, probably not very rugged. Still, to be on the safe 
side, we use the Exchange Monte Carlo method (fTTI) capable 
of, even for glassy systems with many local energy minima, 
sampling the configuration space evenly in a limited time. In 
Exchange Monte Carlo an even number of systems are simu- 
lated in parallel for a sequence of temperatures. After some 
interval, with a probability dependent on the current config- 
uration two systems at adjacent temperatures are exchanged, 
so that one system moves up in temperature, the other one 
down. In our simulations we test for an exchange every 10^ 'th 
sweep. We measure network quantities with the same fre- 
quency as the tests for exchanges. Before we start measuring, 
we run 10^ sweeps for the system to reach equilibrium (which 
is roughly 10 times longer than it takes for all quantities for 
all sizes and temperatures to converge). We use 100 measure- 
ments to calculate intermediate averages. This procedure is 
then repeated for 100 random initial conditions and the aver- 
ages and standard errors of the intermediate averages are the 
values we present below. 

Unless otherwise stated we will use M = 2N. We chose an 
exponential set of 20 temperatures per system size selected in 
a preliminary study to capture the most interesting region. 



III. NUMERICAL RESULTS 

A. Low-energy configurations 

What is the ground state configuration of the YX model? 
First consider a simpler example, — suppose is a multiple 
of three, N = M and the angles are evenly spread out 

over the circle, i.e. 6i = 2m/N. The minimal distance between 



any pair of Oi is the distance In/N between two adjacent ver- 
tices along the circle. There are exact M such pairs, so the 
ground state is E = {(1, 2), (2, 3), • • • , (A^ - 1, A^), (A^, 1)}. Let 
/ = M -h // be the frustration — how far from the lowest 
energy of a model where {6i}iev is unconstrained, the con- 
figuration is. For the state of circularly coupled edges, we 
have (by Taylor expansion, cos0 ^ 1-0^/2 for small 0) 
f = M- Mcos(2n/N) ^ 2Mn^/N^ = 27i^/N for large A^. An- 
other low-energy configuration would be to couple nearby ver- 
tices into A^/3 triangles — E = {(1, 2), (2, 3), (3, 1), • • • , (A^ - 
2, A^ - 1), (A^ - 1, A^), (A^, A^ - 2)}. The frustration in this case 
is twice as large as the circular configuration. If we change 
the system above, so Oi = 27i/N -\- do (for some small pertur- 
bation {) < de < 71 IN) and the rest is the same, then the frus- 
tration increase with 6^ for the circular configuration, while 
it decreases by 67:^^6' IN - for the configuration of isolated 
triangles. This example can be fairly straightforwardly gen- 
eralized to higher edge densities. It suggests that fragmented 
configurations benefit from an irregular distribution of angles, 
while circular distributions should give the lowest energies at 
evenly distributed angles. Since we sample {^J/gv by uniform 
randomness, in the N ^ 00 limit the ground state should be a 
circulant (a graph where vertices are connected to their nearest 
neighbors on a circle). The fragmented states can utilize the 
gaps in the distribution of drawn ^/-values by omitting edges 
across such gaps. For finite sizes, it could happen that the 
ground state is fragmented rather than a circulant. Indeed, 
this is what we see in our simulations. 

In Fig. [T] we show three low-energy configurations for a 
small system size. The lowest energy state, with / = 0.48, is 
fragmented; the other two, with / = 4.14 and 7.95 are more 
elongated, closer to circulants. As we will see, for finite sizes, 
fragmented states like Fig.[TJa) have the lowest energies while 
the circulants have larger entropies making configurations like 
Figs. [TJb) and (c) dominant at higher temperatures. In the 
large temperature limit, the networks are Erdos-Renyi ran- 
dom graphs (EE M)- 



B. Cluster sizes 

Now we turn to the quantitative results. In Fig. [2] we plot the 
relative fraction of the vertices not a part in the largest cluster, 
1 -^1 . Compared with spin models like the XY model, the tem- 
peratures where the interesting behavior occurs is much lower. 
The interesting behavior — the fragmentation of the network 
— is best monitored in a logarithmic temperature scale. The 
onset of fragmentation approaches (as expected from the dis- 
cussion above), zero as the system size increases. In Fig.|2jb) 
we show that for a =1.61+ 0.05 and = 0.22 + 0.03 we have 

1 - ^1 ^ N-^F+iTN''), (3) 

for a smooth function F, and TN^-^^ > 80. The onset of frag- 
mentation thus moves to zero like . Note that this scal- 
ing relation is itself not indicative of critical behavior in the 
rescaled temperature TN^. The observation that the scaling 
region is larger for the smaller than the larger system sizes. 
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FIG. 2 The fraction of vertices not part of the largest connected cluster. In (a) we show the raw data as a function of temperature for four 
different system sizes. In (b) we determine the scaling exponents for low temperatures to a = 1.61 + 0.05 and/? = 0.22 + 0.03. Standard errors 
are smaller than the symbol size and omitted for clarity. Lines are guides for the eyes. 




FIG. 3 The fraction of the second largest cluster, (a) displays the 
raw data of S2 as a function of temperature. In (b) we determine the 
scaling factor y - 1.44 + 0.02. Standard errors are smaller than the 
symbol size and omitted for clarity. Lines are guides for the eyes. 
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FIG. 4 Scaling of the diameter of the largest connected cluster D as 
a function of temperature (a), and the determination of D's scaling 
parameters 6 and e (b). We find 6 ^ 1.52 + 0.02 and e = -0.74 + 
0.02. Standard errors are smaller than the symbol size and omitted 
for clarity. Lines are guides for the eyes. 



which is also not supporting an emergent critical behavior in 

An unexpected phenomenon in the XY model is that there 
seems to be no universality in this size scaling of the network- 
structural quantities. An example of this is seen in Fig. |3] 
where we plot the temperature dependence of the relative 
size of the second largest component ^2- In many systems 
with a fragmentation phase transition, like percolation or net- 
work models of segregation (8; 10), ^2 or S2I si can be used 
to characterize the critical behavior. Also in our case, ^2 
gives a strong signal of the low-temperature fragmentation. 
The heights of the peaks are quite independent of system 
size. The location of the peak scales to zero like T ~ , 
y = 1.44 ± 0.02. This means that the peak of ^2 goes to zero 
slower than the onset of fragmentation as seen in .s-i. 

C. Diameter 

The distance between / and j is the number of edges in the 
shortest path between the two nodes. The largest distance in 
a connected subgraph is the diameter D of that subgraph. If 
the picture from Fig. [T] holds — that at intermediate temper- 
atures are dominated by stretched out configurations, and the 
networks, at lower temperatures, are fragmented — then the 



diameter of the largest connected component should have a 
peak at intermediate temperatures. In Fig. Qa) we plot the 
diameter as a function of temperature. Indeed there is a peak 
moving to zero like T ~ , S = 1.52 ± 0.03, and increasing 
in size as A^"^ with e = -0.74 ± 0.02. For the sizes we test, the 
peak in D occurs at slightly higher temperatures than the peak 
in S2. It also scales slightly faster than the S2 peak towards 
zero. The fact that the peak of D scales sublinearly reflects 
that the intermediate region is a mix of configurations, not all 
stretched out like Fig. [TJb) and (c). This is of course a fun- 
damental aspect of Hamiltonian models — all configurations 
have a finite chance of appearing at all temperatures, but their 
frequencies vary with the temperature. The fact that the peak 
value of D diverges with supports the picture of circulant 
ground states (in the N 00). The peak in Fig.|4jb) becomes 
sharper with larger system sizes, meaning that the increase as 
TN^ is lowered gets more dramatic. This observation is in 
concordance with sl T = phase transition. Our data does, 
however, not give a very strong support for a critical scaling 
of D at the sizes we investigate. 
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FIG. 5 Figures corresponding to Fig. [4jb) for 3N = 2M (a) and 
N - M (b). The scaling exponents for (a) are determined io 6 - 
1.17 + 0.10 and e = -0.52 + 0.02; for (b) they are 6 = 1.25 + 0.10 
and e - -0.61 + 0.02. Standard errors are smaller than the symbol 
size and omitted for clarity. Lines are guides for the eyes. 
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FIG. 6 The number of fragments as a function of temperature (a). 
The solid line is proportional to T'^^^. (b) shows the determination of 
the scaling exponents ^ = 1.94 + 0.07 and 77 = 0.96 + 0.03. Standard 
errors are smaller than the symbol size and omitted for clarity. 



E. Number of isolated subgraphs 

Another feature characterizing the fragmentation is the 
number of connected subgraphs g. In Fig. [6ja) we plot this 
quantity as a function of temperature. The fragmented and 
unfragmented states are rather conspicuously separated and in 
the low-temperature region g = TT'^^^ with the same constant 
r for all system sizes. From Fig.lStb) we see that 



(4) 



g = N-'^G(TN^l 



where f = 1.94 ± 0.07, 7/ = 0.96 ± 0.03, and G is the func- 
tion giving the shape seen in Fig. [6jb). In traditional statis- 
tical physics, for exactly solvable models, the scaling expo- 
nents often have rational values. The g-scaling is consistent 
with a simple picture: if TN^ < 800 we have g ~ T~^^^; if 
TN^ > 1600 we have g - N. To explain the value of we 
note that at the energy close to the ground state is governed 
by angular differences scaling like 1 /N. Hence A^"^ (by the 
leading, square term of the frustration) gives a local energy 
scale, which in analogy to van der Waals' law of correspond- 
ing states implies that T/N~^ is the fundamental quantity for 
the fragmentation, i.e. that ^ = 2. The large TN^ scaling, 
g - N, can be understood from the Erdos-Renyi model that 
is the high temperature limit of the YX model and has the 
same scaling behavior. The low TN^ scaling, g ~ T~^^^, is 
related to the state with fragmented dense clusters. Assume 
that all clusters are of similar sizes, then there are N/g ver- 
tices in each cluster. The number of edges in a dense cluster 
scales like (N/g)^, and thus the total number of edges like 
M = g(N/g)^ = N^/g. Since we assume sparse graphs, 
M ^ N, the number of clusters in such a graph scales like 
g ^ N. By the same argument as above, that TN^ is a funda- 
mental quantity, we get g - T~^^^. 



F. Correlations and radial structure 



D. Density dependence 

We perform most of our analysis for M/N = 2 but will 
briefly touch on how the scaling depends on the density of 
edges. As our example we choose the scaling analysis of 
the diameter seen in Fig. |4jb), but all the other scaling plots 
give the same conclusion, see Fig. [5] Both the networks with 
M/N = 2/3 (a) and M/N = 1 (b) have peaks of the diam- 
eter scaling to zero with increasing sizes. The peak is less 
marked in the sparser graphs in (a) which is natural since they 
are closer to the fragmentation threshold at high temperatures 
(M/N = 1/2). The exponents are also A^-dependent, with S de- 
creasing with decreasing edge density, and e increasing (from 
more to less negative values) as the density becomes larger. 
If the density is larger, e.g. M/N = 4, the scaling we plot, at 
least for the sizes we measure break down (so that there is no 
S such that D peaks at the same TN^ and there is no e such 
that the maximal TN^ is the same for all N). 



To get a more detailed statistical description of the config- 
urations at difl'erent temperatures, we investigate the expecta- 
tion values p of the number of vertices at a distance r from 
a vertex, in Fig. [7] Panels (a) and (b) show the curves be- 
low, and at the maximum diameter. For these curves, the ra- 
dial density decreases exponentially, which, we believe, only 
a fragmentation of the networks can explain. In the T ^ 00 
limit, the radial density curve looks peaked, as expected in 
random graphs (12 ), and similar to observations in large con- 
nected networks (9j. 

In the XY model on two-dimensional lattices, one of the 
central observations is that the spin correlations decay alge- 
braically in the low temperature phase as opposed to an expo- 
nential decay for high temperatures. In Fig. [8j we graph the 
correlation function 

C(r) = |(e'«'''-«^))J, (5) 

where ( • )r denote an average over vertex pairs separated by 
a graph distance r. If all vertex pairs at a certain distance have 
the same relative angular difl'erence C will be one. Comparing 
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FIG. 7 The average density of vertices at a distance r from a vertex for three different temperatures (indicated in the panels) and N = 1600. 
To facilitate comparison, we let the panels cover the same ranges. The widths of the lines are at least the standard error. The inset in panel (c) 
is a blow-up. 
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FIG. 8 The correlation function as a function of r from for three different temperatures and N = 1600, plotted in the same way as Fig. [7] 



the four panels we note that (close to r = 0) the correlations 
decay slowest for the lowest temperatures. This is the same 
behavior as in any classical spin model of statistical physics. 
But in panels (b) and (c) there is a second peak. We understand 
this by analogy to the Erdos-Renyi model — if 2M > N most 
vertices of this model is connected into a "giant component", 
while the rest of the graph consists of small components. In 
our case, there are also disconnected subgraphs for larger tem- 
peratures (i.e., when the graph has a giant component). The 
non-giant components are sparser than the giant and should 
therefore be more volatile in structure. This leads to lower 
correlations at distances close to the diameter of the smaller 
components. 

IV. SUMMARY 

We have investigated the YX model — the XY model where 
the connections, not the spins, are updated — and found a 
complex pattern of size scaling. There are, for the sizes we 
study, three regimes: one regime of fragmented configurations 
dominating at low temperatures, an intermediate-temperature 
regime dominated by stretched-out networks, and a high tem- 
perature, disordered region. The different quantities we in- 
vestigate, all capturing different aspects of the network struc- 
ture, scale to zero with different exponents. Thus the borders 
between the regimes also scale to zero temperature. For ex- 
ample, we measure the diameter (quantifying how elongated 
the largest connected component is) and the size of the largest 
and second largest components. The peak in diameter scales 



to r = like ~ A^-i-52±o.02 ^^iHq the peak of the size of the 
second largest component scale to T = like ~ /^-i -44*0.02 
This multiscaling implies that the picture of three different re- 
gions will change as the sizes increase beyond the ones we 
sample — e.g., the peak of the diameter will, for very large 
sizes, not coincide with the peak of the size of the second 
largest component. 

We do not derive the YX model as a model of some spe- 
cific system in nature or society, but can we a posteriori iden- 
tify such systems? Closest, we believe, is to interpret the YX 
model as describing a social system where the spins repre- 
sent opinions. As we have described the model, it does not 
describe the time evolution of the system. As most interest- 
ing questions about social systems regard their evolution, their 
models need to be dynamic. Dynamics can easily be added, 
most easily by taking the time evolution of e.g. Metropolis 
Monte Carlo sampling ( 21 ). Whether or not this is reasonable 
depends on the system to be modeled. A condition on the sub- 
ject system is that the social ties evolve faster than the opin- 
ions (as in models of social segregation (23) or in some limits 
of coevolution models of networks and opinions dSlfTOl)). 

In the YX model, the role of the N ^ oo limit is, judging 
from our observations, that regime of fragmented states dis- 
appear. There are apparently no emergent singularities. On 
the other hand, the finite-size scaling shows a complex be- 
havior with different scaling exponents. We cannot rule out 
a scenario our scaling parameters converge sls N ^ oo and 
there is a unique parameter;^ such that all quantities signal the 
elongated-configuration regime at a critical TN^. This, we be- 
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lieve, would be the most likely scenario of a phase transition. 
But nothing in our results suggests this would happen. On the 
other hand, the size scaling itself is highly complex. Viewed 
in this way, we have a statistical physics model where the fi- 
nite sizes are more interesting than the thermodynamic limit. 
Moreover, since many real systems (especially in interdisci- 
plinary physics) have restricted sizes, we believe focusing on 
size scaling rather than extrapolating to infinite sizes is a fruit- 
ful future direction for the analysis of statistical-mechanics 
models. 
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